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ABSTRACT 

Recent N-body simulations have shown that the assembly history of galac- 
tic halos depend on the density of large-scale environment. It implies that the 
galaxy properties like age and size of bulge may also vary with the surround- 
ing large-scale structures, which are characterized by the tidal shear as well as 
the density. By using a sample of 15, 882 well-resolved nearby galaxies from the 
TuUy Catalog and the real space tidal field reconstructed from the 2Mass Red- 
shift Survey (2MRS), we investigate the dependence of galaxy morphological type 
on the shear of large-scale environment where the galaxies are embedded. We 
first calculate the large scale dimensionless overdensities (5) and the large-scale 
ellipticities (e) of the regions where the Tully galaxies are located and classify the 
TuUy galaxies according to their morphological types and create subsamples se- 
lected at similar value of 5 but span different ranges in e. We calculate the mean 
ellipticity, (e), averaged over each subsample and find a signal of variation of (e) 
with galaxy morphological type: For the case of 0.5 < 6 < 1.0, the ellipticals are 
found to be preferentially located in the regions with low ellipticity. For the case 
of —0.3 < 5 < 0.1, the latest-type spirals are found to be preferentially located 
in the regions with high ellipticity. The null hypothesis that the mean elliptic- 
ities of the regions where the ellipticals and the latest type spirals are located 
are same as the global mean ellipticity averaged over all types is rejected at 3cr 
level when —0.3 < 5 < 0.1. Yet, no signal of galaxy-shear correlation is found in 
the highly overdense/underdense regions. The observed trend suggests that the 
formation epochs of galactic halos might be a function not only of halo mass and 
large-scale density but also of large-scale shear. Since the statistical significance 
of the overall trend is low, it will require a sample of at least 100, 000 galaxies to 
verify the existence of this correlation. 



Subject headings: cosmology: theory — galaxies: observation 
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INTRODUCTION 



It has been known for long that the physical properties of galaxies such as morpho- 



ter and so on are functions of their environments (Dresslerl 


1980|: 


Postman & Gellei 


1984 




Whit 


more et al. IQQsl: Lewis et al. 20021: 


Gomez et al. 20031: 


Goto et al. 2003: 


Roias et al 


2005 


Kuehn & Rvdedboosl: Blaton et al. 


20051: Bernardi et al 


. 20061; 


Choi et al. 


2007:1 


^ark et al. 


2007 


). Most previous works have largely focused on the galaxies located in the highly dense 



regions and usually quantified the environmental dependence of galaxy properties in terms 
of cross-correlations with local density on small scale, which is believed to be established by 
environment-dependent processes like galaxy-galaxy interaction. 

The currently popular ACDM model predicts that the galaxy properties are correlated 
with not only the small scale but also the large-scale environment. Here the small scale 
and the large-scale environments represent the surrounding regions smaller and larger than 
5/i~^Mpc, respectively. Recent high-resolution N-body simulations of a ACDM cosmology 
have demonstrated that the assembly history of g alactic halos of a given mass change with the 
density field smoothed on sufficiently large scale (IGao et al.ll2005r). Since the galaxy conten t 
depend strongly on assembly history of its host halo (ISpringel et al.ll2005l : ICroton et al.ll2007l ) , 
this numerical finding suggests that the intrinsic properties of galaxies may be correlated 
with the density of large-scale environments. 

In fact, several groups have already found observational evidences fo r the existence 
of cro ss-correlations between galaxy properties and large-scale tidal shear. iNavarro et. al. 



( 120041 ) found that the spin axes of spiral galaxies near th e Local Supercluster lie preferen 



tially on the Supergalactic Plane. iTrujillo et. al.l (120061 ) also showed by analyzing recent 



observational data from large galaxy surveys that the spin axes of void galaxies tend to 
be inclined to the voi d surf aces, keeping the initial memory of tidally induce alignments. 
Pandey fc Bharadwajl (120061 ) have shown that the galaxy luminosity and colour are depen 



dent on the filamentarity of large-scale environment. More recently, iPandey &: Bharadwaj 



( I2OO8I ) have also presented that the star formation rate of galaxies are cross- c orrela ted with 
the filamentarity of the surrounding large-scale structures. [Hernandez et al.l (120071 ) found a 
one-to-one correspo ndence between the value of the gala xy spi n parameter and the morpho- 
l ogica l type, while [Hernandez &: Cervantes-Sodi I (120061 ) and ICervantes-Sodi fc Hernandez 
( 120081 ) have noted that the spin parameter of the observed galaxies is strongly correlated 
with their color and chemical abundance. 

The existence of cross-correlations between galaxy properties and density of large-scale 
environment implies that the galaxy content still have the memory of the initial conditions 
of the Lagrangian regions where the host halos formed. In the scenario of the biased galaxy 
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formation, the galaxy site s correspond to the high peaks of the hnear density field (IKaiser 



1984J : iBardeen et al.lll986l ). But, since the initial density peaks are not spherical, they should 
be characterized not only by the peak heights but also by their ellipticity and prolateness. 
Hence, if the galaxy properties are still linked to the large-scale density, then it is likely that 
they are also linked to the large-scale ellipticity and prolaticity. 

Our goal here is to test this statistical link. Instead of the ellipticity of the large-scale 
density field, however, we consider the ellipticity of the large-scale potential field (i.e., the 
tidal shear). Since the potential field is smoother than the density field, we believe that the 
large-scale potential field refiects more directly the linear conditions. Furthermore, the tidal 
shear field, defined as the second derivative of the potential field, has a large-scal e coherence. 



resulting in highly anisotropic spatial distribution of galaxies at present epoch (IBond et al. 



19961 ). Therefore, the tidal shear field is more readily measurable from the large-scale spatial 



distribution of present galaxies. 

In previous approaches, however, the correlation between galaxy properties and the 
shear of large-scale structure was measured in redshift space. Since the redshift distortion 
effect could contaminate the measurements of the filamentarity of the large scale structure, it 
will be desirable to measure it in real space. Moreover, true as it is that in the linear regime 
the density and the tidal shear at a given region are mutually independent, it is likely that 
the two quantities have developed cross-correlations in the subsequent non-linear regime. 
Therefore, to find an independent relationship between galaxy properties and large-scale 
tidal shear, one has to first account for the density-shear correlations and remove its effect. 

We attempt here to measure observationally an independent relationship between the 
morphological types of nearby large galaxies and the real-space tidal shear reconstructed from 
all sky survey. The organization of this paper is as follows: In §2, the observational data are 
described. In §3, the shear of environment is defined and its mean value averaged over each 
galaxy subsample at similar density is measured as a function of galaxy morphological type. 
In §4, the results are discussed and a final conclusion is drawn. 



2. OBSERVATIONAL DATA: AN OVERVIEW 



A total of 35, 000 galaxies with spectroscopic information a re listed in the TuUy Catalog , 
which was constructed from the ESO/Uppsala full-sky survey (jNilsonlll974j : lLaubertslll982l ). 
Only 15, 922 galaxies in the TuUy catalog have detailed information on their morphological 
types, which were determined by B. Tully according to the conventional Hubble classification 
scheme. For a detailed description of the Hubble classification scheme, see Table 2 in the 
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Third References Catalog of Bright Galaxies (RC3 Ide Vaucouleurs et al.lll99ll ). As for the 
rest 19, 078 galaxies in the TuUy catalog, the stages of their morphological types are described 
as either "irregulars" or "uncertain" or "doubtful" or "spindle", "outer ring" or "pseudo- 
outer R", or "irregulars". These are thus excluded from our analysis. 

According to the morphological type described in RC3, we divide the selected Tully 
galaxies into six samples: E, L, SI, SII, SIII, and SIV. Table[T]lists the class, the morphological 
type (as in RC3), and the number of the galaxies (Ng) belonging to the six samples: The 
sample "E" contains only ellipticals with the Hubble type of cE,EO, EO-1, E+; The sample L 
contains only lenticulars with the Hubble type of SO~^,SO'^ ,S0^; The sample SI consists of 
the spirals of types SO /a and Sa which have tightly wound arms and a largest bulge. While 
the spirals of types Scd, Sd, Sdm and Sm belong to the sample SIV , which have no well 
developed bulge and loosely wound arms. The sample SII contains the spirals of type Sa 
as well as the spirals of type Sab that are intermediate between Sa and Sh. Likewise, the 
sample SIII contains the spirals type of Sh and the spirals of type She. As can be noted, 
the size of bulge and the tightness of the spiral arms decrease from SI to SIV . 

The real space density field was originally r econstructed on 64^ grids in a box of linear 
size 400/i^^Mpc from the 2Mass Redshift Survey JErdogdu et al.lbood . 2MRS '). Basically, it is 
the lin ear density field smoothed with a wiener filter assuming a linear bias. iLee fc Erdogdu 
( 120071 ) used the 2MRS density field to recons t ruct t he tidal shear field on the same 64^ 
grids. As described in detail in lErdogdu et al.l (120061 ) . the real space positions were recov- 
ered by applying the Wiener reconstruction algorithm to the 2MRS data. This algorithm 
basically deconvolve the linear redshift distortions in the radial direction using a distortion 
matrix. For a detaile d description of the method used to recover the real-space position, see 
Erdogdu et aP J2006h . 



The reconstructed tidal shear field consists of a set of three eigenvalues, {Ai, A2, A3}, of 
the tidal shear tensor assigned to each grid. Figure [1] plots the probability density distribu- 
tions of the three eigenvalues of the 2MRS tidal field, assuming Ai > A2 > A3. As can be 
seen, the distribution of the second largest eigenvalue, p{\2)i is almost symmetric around 
the peak at zero with narrow width, while the other two distributions p(Ai) and p{\z) are 
asymmetric around zero with broader shapes. 

The differences among the eigenvalues are related to the shear of environment at a given 
grid. If Ai = A2 = A3, it means that the shear of environment is zero. The larger the mutual 
differences among the eigenvalues are, the higher the shear of environment is. Figure [2] 
plots the 2MRS tidal shear eigenvalue field evaluated on a thin shell at 80/;,~^Mpc shown in 
Supergalactic Aitoff projection. The top, middle and bottom panel shows the Ai, A2 and A3 
field, respectively. In each panel, the darker regions correspond to the low eigenvalues while 
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the bright regions to the high eigenvalues. 



ation method (IHockney &: Eastwood! Il988l ) 



By applying the cloud-in- cell (CIC) interp o ^ 
to the 2MRS tidal shear field, iLee &: Erdogdul (120071 ) also calculated the tidal tensors and 
their eigenvalues at the positions of the selected 15, 922 Tully galaxies. Using the recon- 
structed 2MRS tidal shear field, we measure the values of the shear at the positions of the 
Tully galaxies. In §3, we will see how this measured shears are related to the morphological 
types of the Tully galaxies. 



3. Galaxy-Shear Correlations 

3.1. Variation with the LSS Type 

Let us consider a position, x, where a selected Tully galaxy is located. It depends on 
the signs of Ai, A2, A3 at x, whether the selected galaxy belongs to a void or a sheet or a 
filament or a halo. If all three eigenvalues are positive at x, then the given galaxy is located 
in a halo-like region; If A2 > 0, A3 < 0, then it is in a filament-like region; If Ai > 0, A2 < 0, 
then it is in a sheet-like region; If Ai < 0, then it is in a void-like region. 

We investigate the dependence of the relative abundance of the sample galaxies located 
in halo-like regions on the morphological type, by measuring the conditional number density 
(CND): 

6N,n ~ Ntype ^ AN^n 
Here, iVtype is the number of the galaxies belonging to a given sample, AA'type is the number 
of those galaxies in the sample which are located in halo-like regions, Ng^\ is the number of 
all selected Tully galaxies and AA'^aii is the number of those Tully galaxies which are located 
in halo- like regions. If CND > 1, then the galaxies in a given sample have stronger tendency 
to be located in halo-like regions than the parent sample. 

The CND of the galaxies located in void-like, sheet-like, filament-like regions are also 
calculated in a similar manner. Figure [3] plots the CNDs as histograms for the six samples 
with Poisson errors. As can be seen, the six histograms show all distinct behaviors: The 
CND of the sample E is highest in halo-like regions and lowest in void-like regions; The 
CND of the sample L is lowest in sheet-like regions; The CND of the sample SI is highest in 
filament-like regions; The CND of the sample SII is almost uniform; The CND of the sample 
SIII is highest in void-like and sheet-like regions; The CND of the sample SIV is highest in 
void-like regions. It is interesting to see that in void-like regions the CND of the lenticulars 
have relatively high value compared with that of the ellipticals. 
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Although the result shown in Fig. [3] shows clearly that the galaxy morphology is cor- 
related with the shear of the surrounding large-scale structure, it should not be directly 
translated into an observational evidence for the existence of an independent relationship 
between galaxy morphology and large-scale shear. Since the signs of the three eigenvalues 
of the tidal shears are correlated with the local density, the observed signals shown in Fig. [3] 
could be resulted from the well known correlations between the local density and the galaxy 
morphology. For instance, the observed tendency that the ellipticals are more abundant in 
the halo-like regions may be due to the fact that the local density is usually higher in the 
halo-like regions since the three eigenvalues are all positive and the ellipticals are preferen- 
tially located in the high-density regions. Therefore, to detect the true relationship between 
galaxy morphology and the shear of large-scale structure, one has to use controlled subsam- 
ples where the correlations between galaxy morphology and density are removed. In the 
following subsection, we pursue this task. 



3.2. Variation with the EUipticity 

The shear of environment at a given galaxy position x is caused by the anisotropy in 
the distribution of the surrounding large scale structure which in turn induces asphericity 
in the gravitational pot ential $(x). The asp hericity of $(x) can be quantified in terms of 



ellipticity, e, defined as (IBardeen et al.lll986l ) 



^ A1-A3 

m - 

Here the sum of the three eigenvalues of the tidal shear equals the dimensionless overdensity, 
6 = {p — p)/p where p is the background density: 5 = Yl^=i ^i- If ^^e three eigenvalues, 
Ai, A2, A3, at x have the same value, the iso-potential surface at x has a spherical shape 
with 6 = 0. The larger the differences between Ai and A3 are, the higher the value of e is. If 
the spatial distribution of galaxies located in a region is highly anisotropic, then the value 
of e of the region will be high. In other words, the region has a high ellipticity. The higher 
the degree of the anisotropy in the spatial distribution of the galaxies in a region increase, 
the higher value of the ellipticity the region has. 

At the position of each Tully galaxy, we determine the values of 6 and e. We first 
calculate the means of A6 = 6 — 6g averaged separately over each sample listed in Table 
[H where 6g is the global mean value of the dimensionless overdensity averaged over all 
selected Tully galaxies. If A6 > 0, there is a tendency of the sample galaxies to be located 
in high-density regions. If Ad < 0, the sample galaxies tend to be located in low-density 
regions. If A6 = for all six samples, there is no correlation between galaxy morphology and 
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large-scale density. Figure H] plots ( A(5) of the six samples. As can be seen, there is indeed a 
strong dependence of galaxy-morp hology on the large-scale density, which is consistent with 



numerical results (iGao et al.ll2005l ). 



Now, we calculate the means of Ae = e — ec averaged separately over each sample, 
under the constraint that the value of 6 is fixed in some narrow range to remove the effect 
of morphology- density correlations. Note that is averaged in the same constrained range 
of S. If Ae > 0, there is a independent tendency of the sample galaxies to be located in 
the high-shear regions. If Ae < 0, the sample galaxies tend to be located in the low-shear 
regions. If Ae = for all six samples, there is no correlation between galaxy morphology 
and large-scale shear. 

In Fig. Owe plot (Ae) of the six samples in the top panel for the case where the value of 
6 is fixed in the range of [—0.3, —0.1]. The errors, a, are calculated as one standard deviation 
in the measurement of the mean value: a = ^[((Ae)2) - {Aey]/{Ng - 1) where Ng is the 
number of the sample galaxies selected at the given density range. To demonstrate that 
the correlation between morphology and density is effectively removed by constraining the 
value of 6 to this narrow range of [—0.3, —0.1], we also plot (A(5) of the six samples in the 
bottom panel. As can be seen, the value of (A5), is almost uniform over the six subsamples, 
indicating that the morphology-density correlations are controlled to a negligible level with 
this constraint of —0.3 < 6 < —0.1. We detect a 2.9a signal that the mean ellipticity of 
the sample SIV deviates from the global mean ellipticity. This result suggests that the 
latest-type spirals tend to be located in the regions of high-ellipticity. 

Fig. [6] plots the same as Fig. [5] but with a different constraint of 0.05 < S < 0.1. A 
3(7 signal of morphology-shear correlation is also detected for the sample SIV galaxies. It 
is consistent with the previous result that the latest-type spirals tend to be preferentially 
located in the high-shear environment. Fig. [7] plots the same as Fig. Obut with S in the 
range of [0.2, 0.4]. We detect a 2.2a signal that the mean ellipticity of the regions where the 
galaxies of the sample E are located is lower than the global mean ellipticity, and a 2.4(T 
signal that the mean ellipticity of the regions where the galaxies of the sample SIV is located 
is higher than the global mean ellipticity. Fig. [S] plots the same as Fig. [S] but with 6 in the 
range of [0.5, 1.06]. There are detected a 3. la signal that the mean ellipticity of the regions 
where the sample E galaxies are located is lower than the global mean ellipticity, and a 2a 
signal of correlation between latest type spirals (SIV) and high-shear region. The results 
shown in Figs. [7] and [8] indicate consistently that the ellipticals tend to be preferentially 
located in the low-shear environment while the latest-type spirals tend to be preferentially 
located in the high-shear environment. 



Nevertheless, it has to be noted that the galaxy-shear correlation signal is found to be 
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statistically significant only in bins containing E and SIV galaxies. Therefore, when the 
results over all six bins are considered the null hypothesis of no galaxy-shear correlation is 
still quite acceptable. In the given density range, the null hypothesis is found to be rejected 
at only the ~ 10% confidence level when the results over all six bins are considered. Since we 
consider only those galaxies at similar densities, the sample size, Ng, is quite small, resulting 
in large statistical errors. It is also worth mentioning that neither in the highly underdense 
region with 6 < —0.5 nor in the highly overdense region with 6 > 1.0, no signal of galaxy- 
shear correlation is found since the degree of density-shear correlations is too high to be 
removed in these regions. 



4. DISCUSSION AND CONCLUSION 

By measuring observationally the mean ellipticities of large-scale structures where the 
nearby galaxies are embedded as a function of galaxy morphological types, we have tenta- 
tively detected an independent relationship between galaxy morphology and shear of envi- 
ronment: In the mildly overdense environment with large scale dimensionless overdensity 
of 0.05 < S < 1.06, the ellipticals are found to prefer the low-shear region; In the mildly 
underdense environment with large scale density of —0.3 < S < —0.1, the latest-type spirals 
(Scd-Sm) are found to prefer the high-shear region. No correlation signal, however, is found 
either in the highly overdense {6 > 1) or in the highly underdense environment {6 < —0.5). 

This observational results may be explained by the dependence of halo formation epochs 
on the shear of large-scale environment. It has been known that the formation epochs of 
galact ic halos depend not only on their mass but also on the density of large-scale environ- 



ment (IGao et al.ll2005l ). Here we argue that there is some observational evidence to support 
the hypothesis that the formation epochs of halos also depend on the shear of the large- 
scale environment. In the high-shear environment the halos formed relatively early without 
growing at late time due to strong tidal disruption from the surrounding matter distribution. 
Thus, the late-type spirals are likely to be located in this high-shear environment. Meanwhile 
in the low-shear environment the halos formed relatively recently since there is no strong 
tidal distribution, having grown at late times through hierarchical merging. Hence, the ellip- 
ticals are likely to be found in this low-shear environment. In the highly overdense regions, 
however, even in case that there is tidal disruption from the surrounding matter distribution, 
the other stronger effects of environmental processes should mask the shear-dependence of 
galaxy properties. Likewise in the highly underdense regions, the tidal field is too weak 
to produce any significant galaxy-shear correlation. It will be interesting to test against 
N-body simulations whether or not the formation epochs of galactic halos with similar mass 
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at similar density is a function of the ellipticity of large-scale dark matter distribution. Our 
future work is in this direction. 

It should be noted here that our results suffer from low statistical significance due to 
the small sample size. We have seen only a marginal effect of the large-scale shear for the 
latest-type spirals (Scd-Sm). Meanwhile no signal has been found for the Sb-Sc galaxies even 
though both types of galaxies occupy similar regions. To remove the effect of the stronger 
galaxy-density correlations, we had to constrain the overdensities of galactic regions, which 
resulted in large errors. The null hypothesis that there is no galaxy-shear correlation over the 
six bins is found to be rejected at only the 10% level due to large errors. Given our result, 
it will require a sample of more than 100, 000 galaxies to test the hypothesis at > 90% 
confidence level. This line of investigation will provide new insight into the formation and 
evolution of galaxies in a filamentary cosmic web. 
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Table 1. 



Subsample 


Class 


lype 




E 


Ellipticals 


cE,EO,EO-l,E+ 


1008 


L 


Lenticulars 


so-,so<',so+ 


2532 


SI 


Spirals 


SOa,Sa 


1761 


SII 


Spirals 


Sab,Sb 


3175 


SIII 


Spirals 


Sbc,Sc 


4507 


SIV 


Spirals 


Scd,Sd,Sdm,Sm 


2679 
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Fig. 1. — The probability density distributions of the three eigenvalues of the 2 MRS tidal 
shear field. 




Fig. 2. — The Aitoff Projections of the tidal shear eigenvalue field in Supergalactic coordi- 
nates at 8000 km/s: Ai (top), A2 (middle), A3 (bottom) 



-16- 



a 



0.5 



0.5 



E 



SII 



SIII. 



0.5 



SI __ 



1-^ 



SIV_ 



Void Sheet Fila Halo Void Sheet Fila Halo 



Fig. 3. — The conditional number densities of the galaxies located in the void-hke, sheet-hke, 
filament-like and halo-like regions for the six subsamples with Poisson errors. 
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Fig. 4. — Mean of the density difference averaged over the regions where the TuUy galaxies 
of each subsample are located. The TuUy galaxies are classified into six subsamples accord- 
ing to the galaxy morphological type. The errors represent one standard deviation in the 
measurement of the mean values. 
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Fig. 5. — Mean of the ellipticity difference (top) and the density difference (bottom) averaged 
over each subsample with S in the range of [—0.3, —0.1]. 
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Fig. 6. — Same as Fig. [5] but with 6 in the range of [0.05, 0.1]. 
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Fig. 7. — Same as Fig. [5] but with S in the range of [0.2, 0.4]. 
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Fig. 8. — Same as Fig. [5]but with 5 in the range of [0.5, 1.06]. 



